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ABSTRACT 



["T I . The giant molecular cloud G216-2.5, also known as Maddalena's cloud or 

00 ' the Maddalena-Thaddeus cloud, is distinguished by an unusual combination of 

high gas mass (1 — 6 x 10^ Mq), low kinetic temperatures (10 K), and the lack 
of bright far infrared emission. Although star formation has been detected in 
neighboring satellite clouds, little evidence for star formation has been found in 
the main body of this cloud. Using a combination of mid-infrared observations 
'"Q_ii with the IRAC and MIPS instruments onboard the Spitzer space telescope, and 

Q I near-IR images taken with the Flamingos camera on the KPNO 2.1-meter, we 

identify a population of 41 young stars with disks and 33 protostars in the center 
ci '. of the cloud. Most of the young stellar objects are coincident with a filamentary 

structure of dense gas detected in CS (2-^1). These observations show that the 
^ , main body of G216 is actively forming stars, although at a low stellar density 

00 ' comparable to that found in the Taurus cloud. 
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Introduction 



The weak CO emission f rom t he cold giant molecular cloud G216-2.5 was first discovered 
by iMaddalena fc Thaddeud (119851 ) in a survey of the O rion-Monocerous region. The c loud 
is unusual because it ha s the typical siz e (250 x 100 pc, Maddalena &: Thaddeus 19851 ) and 



mass d - 6 X 10^ M 



0) 



Lee et al.lll994l ) of a giant molecular clou d (hereafter: GMC), yet 
it has a kinetic temperature of only 10 K ( Williams fc Blit3ll998l ). The low ter nperature is 



unusu al for GMCs and is more typical of dense cores in the Taurus dark cloud (jjijina et al. 



19991). Perhaps the mos t unus ual characteristic of the cloud is the lack of star formation. 
Maddalena fc Thaddeud (119851 ) found no clear evidence for star formation within the G216- 
2.5 GMC (hereafter: G216); a result supported by subsequent IRAS maps show ing a distinct 
absence of emission from dust heated by internal young stars (jhee et al.l 1996 ). Considering 
the current distance estimate of 2.2 kpc (iLee et al.lll99ll ). the absence of bright far-IR emis- 
sion only rules out the presence of young high mass stars. Young low to intermediate mass 
stars could have escaped the detection by IRAS. Nevertheless, the evident lack of massive 
star f ormation is in itself unus ual since GMCs almost ubiquitously contain young massive 
stars JwiUiams fc McKeel ^9971 ) . 



Because of the lack of evident star formation, G216 has been considered the best example 
of a quies cent GMC and potentially, a rare example of a GMC before t he onset of star 
formation (IMaddalena fc Thaddeual 19851 ). In the 1990s, near-IR imaging by lLee et al.l (119961 ) 
revealed two small groups of young stars in satellite clouds on the northern edge of G216 
(Fig. [1]). One of these regions is associa ted with a radio source/HII region which was first 
noted by IMaddalena &: Thaddeud (Il985l ). but Maddalena et al. could not rehably associate 
the radio source with the satellite cloud. Near-IR observations of the main body of the cloud 
failed to detect any stellar groups similar to those in the satellite clouds. 



Lee et al.l (119941 ) argued that the main body of the cloud was the site of past star 



formation. Their evidence was the the presence of shell-like structures in the cloud, the 
existence of a large scale velocity gradient across the cloud, and a virial mass five times 
the mass estimated from the luminosity of the ^^CO emission; all of these features could 
be the result of previous episodes of st ar formation which stirr e d up motions in the cloud. 
In support of this idea, HI mapping by lWilliams &: Maddalena! (|l996l ) showed that G216 is 
linked by a 50 pc photodissociation region to two neighboring molecular clouds containing 
young OB stars. In total, these observations have established that G216 is part of a larger 
complex; however, it is not clear whether G216 is a relic of past star formation or a quiescent 
component of a larger cloud complex containing star forming clouds. To date, the relic 
scenario has one major drawback: the lack of a detectable population of stars in the cloud. 
Interestingly, the only sites of previously known star formation directly associated with G216 
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are found in satellite clouds near the surface of G216 that faces the Sh287 star forming region 
(Fig. [1]). This geometry suggests that star formation may have been initiated by interactions 
with Sh287. 

Despite the controversy over the nature of G216, it remains an unusually cold and 
quiescent cloud compared to other GMCs with similar masses. For this reason, detailed 
studies of this cloud may offer new insights into the condition necessary for star formation and 
the factors that drive the rate and efficiency of star formation. We have recently completed 
a Spitzer space telescope survey of the cor e of G216. This core region show s an extended 
region of relatively high column density gas (ILee et al.lll99lUHeyer et al.ll2006l ). but previous 
to this study showed little evidence for ongoing star formation. Using Spitzer and data from 
the KPNO 2.1 meter and SEST telescopes, we show that the central core contains a complex 
of dense cores undergoing star formation. A total of 33 protostars is found, indicating that 
the region is currently undergoing vigorous star formation. In addition 41 young pre-main 
sequence stars with disks are identified. We compare the rate of star formation in G216 to 
that in nearby molecular clouds in Orion and Taurus. 



2. Observations & Data Reduction 



Fazio et al. 


2004. 


Rieke et al. 


2004) 



was mapped with the IRAC instrument in two separate observations on March 25, 2005. 
During each observation, images were taken at series of pointing positions defined by a 5 x 7 
rectangular grid with a spacing between adjancent grid points of 260". This resulted in 
four overlapping but offset maps in the separate IRAC bands providing complete four band 
coverage over a 22' x 31' region. Three dithered 30 second frames were taken at every grid 
position. The high dynamic range mode was selected; each 30 second frame was accompanied 
by a 1.2 second frame. The total integration time for the two observations and three dithers 
was 160.8 seconds. The MIPS data were obtained in four separate observations: two on 
October 10, 2004, one on April 01, 2005 and one on April 08, 2005. In total, four scan maps 
centered on the same position were obtained; these were executed with the medium scan 
rate. Since the direction of the scanning depended on the date observed, the four maps have 
three different scan orientations. This improves the redundancy of the data and minimizes 
artifacts. The scan maps are 25' x 45' in extent with 148" spacings between the scan legs. 

Photometric magnitudes in the four IRAC bands w 'ere extracted with Pho t Vis version 
1.10, an IDL GUI-based photometry visualization tool (iGutermuth et al.ll2004l ). The aper- 
ture size was 2 pixels (2.4") and the sky annulus extended from 2 (2.4") to 6 (7.2") pixels. 
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To convert the signal (measured in DN s ^) to magnitudes, we used zero points of 19.6642, 
18.9276, 16.8468 and 17.3909 in the 3.6, 4.5, 5.8 and 8 fim bands respectively (these are de- 



rived from the calibration of iReach et al.ll2005l ). The signal was multiplied by an aperture 



correction to take into the small aperture size; the adopted values were 1.213, 1.234, 1.379, 
and 1.584 for 3.6, 4.5, 5.8 and 8 /im bands, respectively. The magnitudes from the 1.2 s 
frames were used if the long frame magnitudes were less than 11, 10, 7, and 8 at 3.6, 4.5, 5.8 
and 8 /zm, respectively; at this point, a comparison of the 30 and 1.2 s frame photometry 
showed that the long frames were becoming non-linear. A total of 11,195 sources were de- 
tected in at least the 3.6 and 4.5 bands and 2508 sources were detected in all four IRAC 
bands. 

The MIPS data were processed using the MIPS instrument team's Data Analysis Tool 



(IGordon et al.ll2005l ). which calibrates the data, corrects distortions, and rejects cosmic rays 
during the co-adding and mosaicking of individual frames. The MIPS 24 yum mosaic created 
from the four scan maps has a pixel size of 1.245" and units of DN s"^ pix^^. MIPS 24 fim 
photometric magnitudes were e xtracted using th e IDL implementation of DAOPHOT in the 



IDL Astronomy User's Library (lLandsmanlll993l ) . Sources were first identified with Phot Vis 



An initial determination of their magnitudes was then obtained by aperture photometry with 
PhotVis; the aperture was set to 5 pixels (6.225") and the sky annulus extended from 12 
(14.94") to 15 (18.675") pixels. The zero point magnitude for the mosaicked 24 /xm data 
was 16.05. This assumed a mosaic pixel that was 1/2 the diameter of an instrument pixel, 
a calibration factor of 6.4 x 10^^ Jy per DN s^^ per inst. pix, a zero magnitude flux of 



7.17 Jy, and a correction from 12 pixels to infinity of 1.146 (lEngelbracht et al.l 120071 ). Since 
we measured our photometry with a 5 pixel aperture, a correction from a 5 pixel to 12 pixel 
aperture was determined from the brighter stars in our sample; the correction of -0.428 mag 
was included in the zero point magnitude. A PSF with a 15 pixel radius was then generated 
using 10 bright, isolated stars scattered around the 24 micron image; these were selected 
to be relatively free of contamination from nearby sources and nebulosity. The aperture 
photometry of these 10 stars set the calibration of the PSF fitting photometry. The PSF 
was fit to all the point sources identified in the image using a 2 pixel fitting radius (i.e. only 
pixels in an inner 2 pixel radius were used in the fit). The PSF fit to each source was then 
subtracted and the above procedure was repeated to identify faint point sources hidden in 
the wings of brighter sources. In the second iteration, 13 additional sources were found. A 
total of 253 sources had detection in all four IRAC bands and the MIPS 24 fim band. 

The 70 /im mosaic created from the four scan maps has a pixel size of 9.9" and units 
of DN s~^. We performed point source identification and photometry on the MIPS 70 fim 
mosaic with Photvis. An aperture size of 16" and a sky annulus extending from 18" to 
39" were used. The adopted zero magnitude was of —1.57; this was based on an aperture 
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correction of 2.07, a calibration factor of 1.6 Jy (DN/s/pix) ^, and a zero flux of 0.778 Jy 



f iGordon et al.l 120071 ). 



Deep near-IR data were obtained with the Flamingos camera on the KPNO 2.1-m 
telescope on November 2, 2004. Flamingos data reduction was performed using custom IDL 
routines, including modules for linearization, fiat-field creation and application, background 
frame creation and subtraction, distortion measurement and correction, and mosaicking. 
Point source detection and synthetic aperture photometry of all point sources were carried 
out using Phot Vis. Aperture photometry was performed using radii of 2" .73, 3" .64, and 6" .06 
for the aperture, inner, and outer sky limits, respectively. The photometry was calibrated 
by measuring the magnitude offset between the instrumental magnitudes derived from the 
Flamingos data and the apparent magnitudes from the 2MASS point source catalog. For the 
i^-band data, the offset was dependent on the pixel location of the point source. A position 
dependent offset was derived by fitting a 3rd order polynomial to the offset as a function of 
row and column number. 

We observed in the CS (2 1) line during 12 nights in January of 1993 with the SEST 
15 m telescope. We used the 3 mm schottky receiver coupled to a 2000 channel acousto-optic 
spectrometer with a 80 kHz resolution. The image sideband response was 26 db below that 
of the signal. The object was observed using position switching, with an integration time of 
60 seconds on source followed by 60 seconds on the reference position. This was repeated 10 
times per position, giving a total on source integration time of 10 minutes per position on 
source. We mapped the cloud using a 40" grid in RA and Dec (1950). The beamsize was 
50". 



The calibration was made using the chopper wheel method (IRohlfs fc WilsonI |2004J ). 
The Orion SiO maser was used for pointing. RMS pointing uncertainties were less than 8". 
The total CS(2 1) flux of the Orion SiO maser was repeatable to within 10%, except on 
the last night. We note that only 21 points were observed on the last night: these points 
were typically non-detections. The average system noise (after correcting for the effect of 
the atmosphere) was 450 K, the maximum was 590 K. The data were reduced using the 
Grenoble CLASS software using standard procedures. A beam effi ciency of 0.75 was used 



to convert the antenna temperatures into main beam temperatures (IHenkel et al.lll994l ): all 
our results are on a main-beam temperature scale. 
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Results 



3.1. Identification and Classification of Young Stellar Objects 

Young stellar objects (YSOs) with infrared excesses due to dusty envelopes and disks 
were identified and classified in G2 16 using methodologies developed in a series of publica- 



tions (for a detailed discussion see lAUen et al.l 120071 . IWinston et al.l 120071 . iGutermuth et al. 



20081 ). Since these methologies were developed for star forming regions in the nearest kilo- 
parsec, they were modified to take into account the greater distance t o G216. In addition, a 
new method for identifying protostars was also implemented (also see lMegeath et al.ll2009l ). 
The overall approach can be decomposed into three steps: the identification of candidate 
YSOs with infrared-excesses, the rejection of extragalactic sources contaminating the sample 
of candidate YSOs, and finally, the classification of the remaining YSOs as either protostars 
or pre-main sequence stars with disks. These steps were implemented using a suite of mid- 
IR color and magnitude criteria. In the following overview of our methodology, we do not 
entirely segregate these three steps since the details depend on the specific colors and mag- 
nitudes employed. 

For the purposes of this paper, we divide our sample of YSOs with infrared excesses 
into two gross evolutionary categories: protostars and young stars with disks. We effectively 
classify objects by the slope of the spectral energy distribution (SEDs). Flat or rising SEDs 
with increasing wavelength (i.e. a > —0.3 where a = dlog\F\/ dlogX) are indicative of an 
infalling envelope; these include Class 0, Class 1 and fiat spectrum sources. Decreasing SEDs 
with increasing wavelengths (a < —0.3 ) are indicative of stars surrounded by circumstellar 
disks. This approach is derived from the pre-Spitzer classification scheme of iGreene et al. 



fll994h 



The starting point of our analysis was the catalog of all point sources detected by 
2MASS, Flamingos, IRAC and/or MIPS; this was compiled from the point source photom- 
etry described in Sec. 2. The first steps were to eliminate point sources which were only 
partially covered in our S'pzteer observations and to reject photometric magnitudes with high 
uncertainties. In this paper, we only consider the 0.16 sq. degree region covered by all four 
IRAC bands as well as the MIPS 24 and 70 /im bands; the Flamingos data cover a somewhat 
smaller region. In the following analysis, uncertainties of < 0.1 mag were required for the 
[3.6] and [4.5] photometry, < 0.15 mag for the [5.8] and [8] photometry, < 0.25 mag for the 
[24] photometry and < 0.15 mag for the Flamingos JHK photometry. Sources were kept in 
the inital point source catalog as long as they contained at least one photometric magnitude 



Using models of protostars and young stars with disks, lAUen et al.l (120041 ) and lMegeath et al 



( I2OO4J ) found the IRAC [3.6]-[4.5] vs. [5.8]-[8] diagram a powerful means for identifying in- 
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frared excess sources (Fig. [2]); these sources include both protostars and young stars with 
disks. We select IR-excess sources which satisfy the following criteria: 

[3.6] - [4.5] > 0.2 + a[3.6]-[4.5], [5.8] - [8.0] > 0.35 + a[5.8]-[8] (1) 

Many sources can be detected in the IRAC 3.6 /xm and 4.5 /xm bands but are not detected in 
the longer wave length IRAC and MIPS band s because of the lower sensitivity to photospheres 
in these bands. iGutermuth et al.l (120041 ) and lWinston et al.l (120071 ) used combined IRAC and 
near-IR color-color diagrams to identify young stellar objects that lack detections longward 
of 4.5 /im. However, given our exceptionally sensitive 5.8 and 8 /xm data (compared to other 
surveys of star forming regions - this results from our relatively long integrations and the 
lack of a bright nebulous background), the comparatively low sensitivity of the 2MASS data 
for this distant region, and the ~ 0.1 mag. scatter in the deeper Flamingos photometry, we 
find the combined IRAC/near-IR diagrams do not show reliable new YSOs. We do find that 
the addition of the 24 /^m data does result in the identification of new YSOs. In particular, 
young stars with disks may not show strong excesses i n the IRAC bands du e to holes in 
their inner disks; these are the transition disks sources (IMuzeroUe et al.l 120041 ). Among the 
sources with detections in the MIPS 24 /im ban d, we search for youn g stars with transition 
disks by using the color criteria (Fig. [21 also see lWinston et al.l 120071 ): 



[8] - [24] > 1 



(2) 



The sample of infrared excess objects sele cted by the above criteria is contaminated by 
galaxies with strong PAH emission and AGN (jStern et al.l 120051 ). PAH dominated galaxies 
exhibit colors distinct from those of YSOs; we eliminate such galaxies u sing the empirically 
derived color criteria given in the appendix of IGutermuth et al.l (120081 ). In contrast, AGN 
show colors similar to those of YSOs and must be distinguished by their fainter magnitudes. 
We determined the optimal threshold magnitude by examining the density of AGN in survey 
data from the SWIRE legacy program; the SWIRE fields were chosen to study extragalacti c 
populations with a minimum of contamination from galactic objects ( iLonsdale et al.ll2003l ). 
We used the point source catalog the from Elias-N2 field obtained from the NASA/IPAC 
Infrared Science Archive. In Fig. [31 we show the number of AGN with colors satisfying our 
YSO color criteria (Equations 1 & 2) as a function of their 4.5 /xm magnitude. The number 
is corrected for the smaller size of the G216 field. Assuming that the density of AGN is 
the same in the SWIRE and G216 fields, the histograms show that the number of expected 
AGN exceeds the number of YSOs for [4.5] > 15. However, instead of ad opting a constant 



magnitude threshold, we use the color dependent magnitude threshold of IGutermuth et al. 
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( 120081 ) ■ We retain the color dependence of the Gutermuth et al. criteria, but we lower the 
magnitude threshold by 0.5 mag to account for the greater distance to G216. Consequently, 
for a source to be classified as an AGN, it must satisfy all of the following criteria: 



[4.5] - [8.0] > 0.5 
[4.5] > 14. + ([4.5] - [8.0] - 2.3)/0.4 
[4.5] > 14 

In addition, a source would have to satisfy one of the following three criteria before it is 
classified as an AGN: 



[4.5] > 14.5 + ([4.5] - [8.0] -0.5) 
[4.5] > 15 
[4.5] > 15 - ([4.5] - [8.0] - 1.2)/0.3 

Sources with colors consistent with Equations 1 and 2 but satisfing the above criteria are 
rejected from the YSO sample. The division between sources selected as YSOs and AGN is 
shown in the [4.5]-[8] vs. [4.5] diagram displayed in Fig.m This figure shows that the adopted 
threshold is typically brighter than 15 mag at 4.5 /im, and consequently, will reject the vast 
majority of contaminating AGN (Fig. [3]). Finally, to minimize galactic contamination in the 
sample of transition disks (which do not have the colors of AGN), we deleted one transition 
disk candidate with [4.5] > 15 mag. 

The above criteria result in a sample of YSOs containing a mixture of protostars and 
young stars with disks. To distinguish between protostars and young stars with disks in 
that sample, and to further increase the sample of protostars, we apply criteria utilizing only 
the three most sensitive bands for detecting protostellar sources: the 3.6 /xm, 4.5 /xm and 
24 ^m bands. Protostars may not be detected in the 5.8 fim and 8 /xm bands because of the 
lower sensitivity in these bands, the bright nebulosity from the strong hydrocarbon features 
between 5 and 9 /im, and the flattening and dip in the protostellar SED apparent in the 
5-10 /im spectral regime. Furthermore, protostars may not be detected in the near-IR bands 
due to their high extinction. 

We identify protostars by adop ting the criteria that objects exhibiting SEDs with spec- 
tral indices a > —0.3 are protostars (jGreene et al. 1994 . Megeath et al. I l2009h . We implement 



this criteria by determining the photometric colors of a power-law SED with a > —0.3 using 
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the IRAC and MIPS spectral response curves posted on the Spitzer Science Center website. 
We find that protostars exhibit the following colors (Fig. [2]): 

[3.6] - [4.5] > 0.652, [4.5] - [24] > 4.761 (3) 

The sample of protostars selected by these color criteria also suffers from contamination 
by galaxies. Specifically, the [4.5]-[24] vs. [24] diagram (Fig. Hj) shows a distinct clump of 
faint sources that satisfy our color criteria; these sources are scattered evenly across the sky 
and are probably extragalactic (Fig. [5]). In Fig. [3l we plot the the number of galaxies in 
the SWIRE field which satisfy the protostar color criteria (Eqn. 3) as a function of 24 /im 
magnitude; once again the number has been scaled to account for the smaller size of the G216 
field. We compare this to the number of protostars candidates identified in the G216 field. 
We find that for [24] < 9, the number of protostar candidates exceeds the number of expected 
galaxies. To minimize contamination from galaxies, we adopt the following empirical criteria 
for YSO: 

[24] < 9.15 (4) 

We find one faint protostar candidate ([24] = 9.1 mag.) clustered with the brighter protostars; 
we have set the threshold for YSOs to 9.15 mag so that this source is included in our YSO 
list. The application Equations 3 & 4 results in a sample of protostars with the following 
three pedigrees: 

1. objects identified as YSOs using the previous critiera but not previously classified as 
protostars 

2. resurrected YSOs which were previousy rejected as AGN but are now considered proto- 
stars (since deeply embedded protostars can have weak 4.5 /xm emission, the rejection 
of galaxies from the protostar sample is better done at 24 /im). 

3. newly identified protostars which were not previously identified as infrared excess 
sources because of a lack of [5.8] and/or [8] photometry. 

Although limited by a lower angular resolution and sensitivity than the other Spitzer 
bands, the 70 fim data adds an additional means to identify cold protostars. One additional 
source is identified as a protostar on the basis of the 70 fim data. This source was not 
initially identified as a protostar since its [3.6]-[4.5] color is similar to a young star with 
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disk. Because of its faintness and its color, it was initially classified as an AGN. However, 
given its brightness in the 70 fim band and its extremely red color, we include this in our 
protostar sample (source 14 in Table 1). The remaining sources with 70 fim band detections 
were not classified as YSOs; these sources are fainter in the 70 fim band than the previously 
identified YSOs and have a 70 /xm magnitude > 3. In comparison, only one previously 
identified YSO has a 70 /im magnitude > 3; this was identified as a protostar through its 
3.6 /zm, 4.5 /xm and 24 /im photometry (source 12 in Table 1). To further test whether the 
unclassified sources with [70] > 3 mag are background galaxies or YSOs, we examine their 
spatial distribution. We find that these faint sources are distributed uniformly in the 70 /im 
mosaic; this is in distinct contrast to the protostellar sources with 70 fim magnitudes < 3 
which are concentrated in the center of the field. We conclude that the unclassified sources 
with [70] > 3 mag are primarily galaxies and do not include them in our YSO sample. 

In total, 33 protostars were identified (Table 1). An altern ative approach to classif y 



protostars using the [4.5]- [5.8] colors of YSOs was proposed by iGutermuth et al.l (120081 ). 
With only three exceptions, the colors of our protostars are consistent with this other criteria 
(Fig. [21). The remaining YSOs with infrared excesses which were not classified as protostars 
were classifed as young stars with disks. A total of 41 stars with disks were identified 
(Table 1). The distribution of protostars and young stars with disks are shown in Fig. [61 

The criteria described above are designed to minimize the contamination of the YSO 
sample by galaxies. The success of these criteria are apparent in Fig. \5[ where we show the 
distribution of the YSOs and the identified contaminating galaxies. While the galaxies are 
uniformly distributed across the field, the YSOs do not exhibit the same uniform distribution. 
Nevertheless, we expect a residual of galaxies in the YSO sample. To determine the level of 
residual extragalactic contamination, we applied the full suite of criteria described above to 
the SWIRE catalog for the Elias-N2 field. After the rejection of likely galaxies, we estimate 
the remaining contribution of extragalactic contamination to our sample is approximately 6.2 
protostars and 7.5 stars with disks. After the subtraction of this background, the corrected 
numbers of bona-fide protostars and young stars with disks are 27 and 33.5, respectively. 
We note that most of the protostars contamination occurs at the fainter magnitudes, with 
half of the extragalactic sources misidentified as protostars having [24] magnitudes between 
8.65 and 9.15 (Fig.[3|). 



3.2. Distribution of Young Stellar Objects and Dense Cores 



In Fig. [6l we show the distribution of young stellar objects overlaid on the IRAC images 
of the region. The MIPS 24 fim image is displayed in Fig. [3 overlaid on the MIPS data 



- 11 - 



are contours of velocity integrated CS (2 1) emission. The extent of the CS map, which 
is smaller than the 24 fim mosaic, is also displayed. The 33 protostars and 41 stars with 
disks are concentrated in a region which extends diagonally through the map. The dense 
gas traced by the CS (2 — > 1) emission forms a clumpy ridge extending 10 pc diagonally 
across the field. Most of the protostars are coincident with detectable CS emission; the six 
exceptions are outside the region mapped with the SEST. Most of the stars with disks are 
also concentrated near the dense gas ridge; the remainder are found to the south of the 
ridge. The region north of the ridge is almost devoid of YSOs. The number of protostars 
outside the ridge is equal to the number of expected contaminating galaxies; this suggests 
that the protostars are confined to the ridge with the sources outside the ridge being likely 
misidentified galaxies. 



The H-K vs. K plot in Fig. H] shows the iBaraffe et al.l (119981 ) pre-main sequence tracks 



for 1 Myr YSOs plotted over the G216 sources. Adopting an age of 1 Myr, most of the pre- 
mains sequence stars would have masses of 0.25-1.2 M0 (Fig. Hj), with the faintest identified 
stars having masses close to the Hydrogen burning limit. A compact red nebula is apparent 
in Fig. El the nebula exhibits bright 8 /im emission indicative of UV heated hydrocarbons. 
The UV source is hkely to be a young star with disk found within the nebula (source 56 
in Table 1). Interestingly, the infrared excess from this star is weak except at 24 /xm. The 
photometry may be affected by the nebular emission; however, visual inspection of the PSF 
subtracted 24 /im mosaic shows that the photometry is not significantly contaminated by 
the nebula. The star is one of the more luminous pre-main sequence stars in this region of 
the cloud. Adopting an age of 1 Myr, the bright source toward the nebula lies well above 
the 1.2 solar mass upper limit for the tracks in Fig. |H Two others star/disks (sources 38 and 
59 in Table 1) have similar near-IR l uminosities. A cornpariso n of the photometry of these 



sources with Herbig Ae/Be stars from lHillenbrand et al.l (119921 ) shows that these sources are 



consistent with A type stars, with K-band magnitudes and H-K colors similar to the A2 star 
HD 245185. Although spectroscopy is needed to confirm spectral types, these data suggest 
that the most massive young stars in the region are early A type stars. This would explain 
why the central core of G216 is not prominent in IRAS maps of the far-IR emission. 



4. Discussion 

Although these observations clearly demonstrate that star formation is ongoing in the 
central core of G216, the 74 YSOs are both small in number and distributed over a region 
of 230 pc^. The high proportion of protostars indicates that this is a young population 



of stars, and is not the result of the past episode of star formation proposed by iLee et al. 
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( 119941 ) ■ However, such an older population may not exhibit the bright infrared excesses used 
to identify the YSOs in this study, and we cannot rule out its existence. 

How does this compare to other molecular clouds? We compare G216 to two of the 
best studied molecular clouds near the Sun: the Taurus dark cloud complex, the prototype 
cloud for distributed star formation, and the Orion A giant molecular cloud, the prototype 
cloud for clustered star formation. Since a comprehensive Spitzer catalog for Taurus has 
not been published, we cannot estimate the number of YSOs with IR-excesses that would 
be selected in a 2200 pc distant Taurus cloud. Instead, we take a recent compilation of all 
known YSOs in Taurus (K. Luhman, P. Com); this contains a total of 296 sources. The entire 
Taurus cloud is extended over a region 1408 sq pc; the YSOs fill only a fraction of this region 
(lUngerechts &: Thaddeud 119871 ). For Orion A, we use the Spitzer Orion su rvey and apply 



simil ar criteria as applied to G216 for identifying YSOs with IR-excesses (IMegeath et al. 



2009h. A total of 19 35 YSOs are identified in the Spitzer survey, which covers 261 sq pc 
jMegeath et~aDl2009h . 



We now compare the densities of YSOs in the G216, Orion and Taurus samples. To 
compare these three samples, we must account for the larger distance of G216. This can 
be done by invoking a cutoff in the J-band magnitude equivalent to the faintest J-band 
magnitude of the G216 YSO sample. From Fig. [HI this value is J — DM < 6.15, where 
DM is the distance modulus. Sources undetected in the J-band, typically prototstars, were 
not eliminated. This is not a rigorous determination of the completeness in the G216 data; 
however, this crude approach is sufficient for our comparison of stellar densities. Using the 
J magnitude cutoff, we estimate that we would detect 1372 sources in the Orion A cloud if 
it were placed at 2200 pc. In Taurus, we would detect 183 sources at 2200 pc. 

In Fig. [HI we show the distribution of local YSO surface densities measured around each 
YSO for G216, Taurus and Orion A. The surface density was calculated using the nearest 
neigh bor method, in which density = n/fTrr^) where r„ is the distance to the nth nearest 
YSO (iGutermuth et al.ll2005l . ICasertano fc Hutlll985l ). This was done for n = 5 and n = 10 
using the total sample of YSOs in each cloud and using the sample of YSOs that satified the 
J — DM criteria. The resulting surface densities in G216 are much lower than the Orion A 
cloud, but are similar to those found in the Taurus dark cloud complex. 

The density and configuration of YSOs in G216 are similar to those in the Taurus 
cloud; with small, low density groups of stars distributed among filamentary clouds. How- 
ever, the cloud mass of G2 16 is ten times tha t of the Taurus cloud and is comparable to 
that of the Orion A cloud (IWilson et al.ll2005l ). Interestingly, the prop erties of the turbu - 
lence in G216 seems to be similar to more active star forming GMCs ( iHeyer et al.ll2006l ). 
It is, however, unlikely that G216 will evolve into an Orion cloud given the estimated star 
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formation rate. Assuming a 0.43 Myr year lifetime for protostars (IHatchell et al.l 120071 ). we 
estimate a star formation rate of 60 stars per Myr for G216. Thus, at this rate, it would take 
17 Myr for G216 to form 1000 sta rs. This is longer than the measured lifetimes of molec- 
ular clouds (IHartmann et al.l 120011 ). G216 could produce a Taurus-like region in 3 million 
years. The CS emission in G216 exhibits lower brightness temperatures and slig htly lower 



linew i dths than those typica lly exhibited by st ar forming dense cores in GMCs (ILee et al. 



1994 IWilliams fc Blitall998l ). iLee et al.1 (119941 ) suggested that the G216 cores are similar 



to those in Taurus; however, a direct comparison of the dense gas is complicated by the low 
spatial resolution (0.6 pc) of the G216 CS (2 1) map. We defer this analysis to a future 
paper using existing multi-transition CS data. 



5. Summary 

We present the first detection of ongoing star formation in the main body of G216, 
otherwise known as Maddalena's cloud or the Maddalena-Thaddeus cloud. We identify 
33 protostars (6.2 of which are expected to be contamination) and 41 young stars with 
disks (7.5 of which are expected to be contamination). This demonstrates that G216 is not 
quiescent, but forming young stars at an approximate rate of 60 stars per million years. The 
most luminous pre-main sequence stars have near-IR luminosities similar to Herbig Ae stars, 
consistent with the lack of massive star formation in the cloud inferred from IRAS far-IR 
maps. The star formation is concentrated in a 10 pc long molecular ridge detected in CS 
(2^1) observations. The density of star formation is much lower than that found in more 
active GMCs such as Orion GMCs, and is similar that in the Taurus dark clouds. This may 
result from the presence of less dense gas and smaller dense cores compared to active star 
forming regions such as Orion. Future work should concentrate on detailed comparisons of 
the molecular cores of G216 with those in Taurus, Orion and other molecular clouds, with 
the goal of understanding how the properties of the gas determines the rate and density of 
star formation. 
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Fig. 1. — CO (1-0) map of G216-2.5 with the lo cation of the Sp itzer map marked in the 
white box; this map w as made from the data of (iLee et aLlll994j ). Areas of star formation 
previously identified by iLee et al.l (Il996l ) are indicated. 



- 16 - 




Fig. 2. — Color-Color Diagrams. All sources with photometry satisfying the uncertainty 
requirements are plotted as dark dots. Those classified as protostars are circles, and those as 
stars with disks are squares . The vectors r e presen t 5 magnitudes of extinction in the iiT-band 
using the reddening law of iFlaherty et al.l (120071 ). 
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Fig. 3. — Top: Histogram of expected AGN contamination as a function of 4.5 //m mag- 
nitude. Solid line: distribution of 4.5 iim magnitudes for our sample of YSOs in G216. 
Dot-dashed line: the estimated distribution for extragalactic sources misidentified as YSOs 
based on the SWIRE data. Bottom: Histogram of expected galaxy contamination to our 
protostar sample. Solid line: distribution of 24 fim magnitudes for our samples of proto- 
stars. Dot-dashed line: the estimated distribution of misidentified extragalactic sources, 
again based on the SWIRE data. 
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Fig. 4. — Color-Magnitude Diagrams. Symbols are same as Fig. 2. The line drawn on the 
[4.5]- [8.0] vs. [4.5] plot represents the criteria used for eliminating AGNs from the data set. 
Note that only stars with disks were required to satisfy this criteria. A large concentration 
of faint, red extragalactic sources is also apparent in the [4.5]-[24] vs. [24] plot (centered on 
[24] ~ 10.5 and [4.5j -[24] ~ 5.5). The H vs. H — K diagram shows the 1 Myr isochrone from 
Baraffe et al. Jl998h : the upper and lower arrows are reddening vectors corresponding to 10 



Ay at masses of 1.2 Mq and 0.25 Mq, respectively. The [3.6]-[4.5] vs [4.5] diagram show the 
photometry of the YSOs relative to all sources detected in the [3.6] and [4.5] bands. 
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Fig. 5. — Distribution of candidate YSOs, AGN, and star-forming galaxies. Green filled 
circles are protostars, green squares are young stars with disks, red asterisks are AGN and 
red Xs are galaxies identified by their their faint 24 fim magnitudes. 




Fig. 6. — Top: the distribution of protostars. Bottom: the distribution of young stars with 
disks. The background image is the combined IRAC 3.6 (blue), 4.5 (green) and 8 /im (red) 
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Fig. 7. — MIPS 24 micron image of the central region containing the young stellar objects. 
The contours show the velocity integrated CS (2 —* 1) emission; the contour levels are 0.75, 
1, 1.25 and 1.5 K km s^^. The extent of the CS map is also shown (delineated by the straight 
lines). The 1 pc scalebar is shown for the adopted distance of 2.2 kpc. 
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Fig. 8. — Histograms of the J-band magnitudes for the Taurus, Orion A and G216 samples. 
The J-band magnitudes have been corrected for the distance modulus (DM), adopting dis- 
tances of 142 pc, 414 pc and 2200 pc for Taurus, Orion and G216 respectively. The faintest 
source in the G216 sample has a J — DM — 6.15 (Table 1). 
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Fig. 9. — Histograms of YSO nearest neighbor densities for G216, Taurus and Orion A. 
The densities are calculated for each YSO using the equation density = n/{7irD where r„ 
is the distance to the nth nearest YSO and n is set to 5 (left) and 10 (right). The upper 
line in each plot outlines the histogram using all the YSOs in each of the three clouds. The 
shaded histograms use only YSOs that satisfy the condition J — DM < 6.15 mag, or are not 



detected in the J-band. The dashed line marks the density of 10 pc ^; iMegeath et al.l (120091 ) 



define groups and clusters as regions with 10 or more YSOs that have a surface density in 
excess of this threshold. Consequently, G216 does not contain a group or cluster by this 
definition. 
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